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Abstract 

A new method, 'Configurationally-Resolved-Super- Transition- Arrays', for calculation of the spec- 
tral absorption coefficient in hot plasmas is presented. In the new method, the spectrum of each 
Super- Transition- Array is evaluated as the Fourier transform of a single Complex Pseudo Parti- 
tion Function, which represents the exact analytical sum of the contributions of all constituting 
unresolved transition arrays sharing the same set of one-electron solutions. Thus, in the new 
method, the spectrum of each Super- Transition- Array is resolved down to the level of the (unre- 
solved) transition arrays. It is shown that the corresponding spectrum, evaluated by the traditional 
Super- Transition-Arrays (STA) method [A. Bar Shalom, J. Oreg, W.H. Goldstein, D. Shvarts and 
A. Zigler, Phys. Rev. A 40, 3183 (1989)], is just the coarse grained Gaussian approximation of 
the Configurationally-Resolved-Super- Transition-Array. A new computer program is presented, 
capable of evaluating the absorption coefficient by both the new configurationally resolved and the 
traditional Gaussian Super- Transition-Arrays methods. A numerical example of gold at tem- 
perature IkeV and density 0.5 gr/cm^ , is presented, demonstrating the simplicity, efficiency and 
accuracy of the new method. 



I. INTRODUCTION 



The radiative opacity is an essential factor governing the structure and evolution of stars 
as well as laboratory plasmas jsl- In plasmas containing medium to high Z elements, 
at least some of the electrons remain bound to the ions even at very high temperatures and 
densities (e.g. iron at the center of the sun). As was first recognized by Edward Tellerj^], in 
part of the spectral range, the opacity in these plasmas is dominated by photo-absorption 
of dipolar radiative transitions between electronic states of the ions (line absorption). 

The present work focuses on line absorption in plasmas in Local Thermodynamic Equilib- 
rium (LTE). Historically, the development of the theory and computational approach to this 
process required a major theoretical effort (e.g. jst, j3]). The atomic states are evaluated, 
in all methods, by a perturbation expansion, using the radial average potential approxima- 
tion as the zero order. In this order, the equation for the many-electron states is reduced to 
equations for the one-electron states in an average radial self-consistent potential due to all 
other electrons. For light elements, Schroedinger equation is sufficient, for heavier elements, 
relativistic treatment is required and the Dirac equation should be solved. The zero order 
many-electron states ("configurations") and energies are characterized by the occupation 
numbers of degenerate groups of one-electron states ("shells"). Mutual electron-electron 
interaction removes the degeneracy and splits the configurational energy. This effect is eval- 
uated as a first order correction, i.e. as the sum of expectation values, in the zero order 
states, of energies due to coulomb interaction between all pairs of electrons. The well known 
Slater integrals represent the radial part of these expectation values. 

The evaluation of the absorption coefficient requires a compromise between the needed 
spectral resolution and the available computer resources (for a representative list of codes 



see e.g. B,^ 



10|). Existing methods may be classified according to the resolution of the 



description of electronic states, and of the contribution of transitions between them to the 
absorption coefficient. The most resolved treatment available is the Detailed Line Account- 



ing (DLA, e.g. 



11|). For complex configurations with many states, the number of transitions 



becomes prohibitive for numerical calculations. In this case, one may turn to the Unresolved- 
Transition-Arrays (UTA) method. In this method the spectral absorption coefficient due 
to a transition array consisting of all single-electron transitions between a specific pair of 
configurations is assumed to be of a Gaussian shape. This method is made practical by 
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the analytical formulae for the three lowest energy-moments (actually cumulants 12|]) of the 



UTA spectrum 



131 ]. in terms of reduced matrix elements of the dipole operator, Slater inte- 



grals and shell occupation numbers. For heavy ions, the UTA method becomes unpractical 
due to the enormous number of possible configurations. The Super- Transition-Array (STA) 
method [l^,[l5|, 16 1 represents a further compromise which allows the evaluation of opacity 
with less computational effort at the cost of spectral coarse graining. The coarse graining 
is obtained by grouping shells, with adjacent energies, into supershells, configurations into 
superconfigurations (SCs) and correspondingly transition arrays into supertransition arrays. 
The relative simplicity of the evaluation of the coarse-grained spectral absorption coefficient 
is based on three assumptions (on top of the UTA assumption): 

a) The basic superconfiguration assumption: All configurations which form a supercon- 
figuration share the same radial potential with the same set of one-particle solutions. 

b) The high-temperature approximation: The spread of the energies of configurations 
within a superconfiguation is much smaller than the plasma temperature. In this limit 
the Boltzmann factor which determines the relative probability for a configuration within 
a superconfiguration may be evaluated to zero order only, i.e. as the sum of single-electron 
energies in the mean potential. Electron-electron interaction energy adds a superconfig- 
urational average factor common to all configurations within a superconfiguration. This 
corresponds to the use of the Gibbs-Bogoliubov-Feynman bound [17] as an estimate for the 
Boltzmann factor. 

c) The unresolved supertransition array assumption: The spectra of all UTAs which form 
a STA merge into a single Gaussian shape. 

With these three approximations the summation of contributions of all UTAs to a STA 



may be performed analytically 



1J],[15|. By the third approximation, one needs only to eval- 



uate the three lowest energy cumulants of the STA spectrum. By the second approximation, 
the relative probabilities of configurations are the same as of those in a system of indepen- 
dent particles in a potential well 18|. This enables the derivation of analytical formulae for 
the moments (and cumulants) in terms of Slater i nteg rals and partition functions jl^ . {igI] 



l|,|20|,|21]. 



which may be evaluated by recursion relations 

In reference [23|], a way to avoid the high-temperature approximation was shown. An 
analytical formula was written for the partition function with the full Hamiltonian including 
electron-electron interaction. This result was enabled by the application of the Hubbard- 
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Stratonovich transformation 



2J| , which ehminates of the quadratic dependence of the energy 



on the shell occupation numbers (at the cost of introduction of an auxiliary random field). 

The approximations in assumptions (a)-(c) are controlled by the choice of the degree of 
spectral coarse graining. In the extreme choice of one shell in a supershell, and one config- 
uration in a superconfiguration, and a different average radial potential with a different set 
of one-electron states separately for each configuration, one reaches the UTA limit. Clearly, 
within the framework of the UTA model, assumptions (a)-(c) are exactly satisfied. The op- 
posite extreme choice is of one supershell consisting of all shells and one superconfiguration 



consisting of all configurations in all degrees of ionization. In reference [IJ] it was shown 
that this choice of maximum spectral coarse graining reproduces the results of the average 
atom (AA) model, presented in reference 25|. 

As mentioned above, for heavy ions, calculations with the UTA resolution are impractical 
due to the enormous number of possible configurations. On the other hand, the AA model 
is too crude since it wipes out the spectral structure observed in experiments . The STA 
method allows for a tune up of the resolution by an iterative refinement procedure starting 
from the AA model, increasing the number of superconfigurations and number of different 
radial potentials with different sets of one-particle solutions. The refinement process is 
stopped when both the values Rosseland and Planck mean free path (MFP)[l] converge to 
constant values. A typical STA run reaches convergence with a few to a few tens of SCF 
solutions per degree of ionization, and many more superconfigurations sharing the same 
potential and set of single particle energies and orbitals. Actually, as will be demonstrated 
by the numerical example in the present work, when a strict convergence test based on the 
spectral details rather than the integrated values of Rosseland and Planck MFP is imposed, 
convergence of the STA refinement process in heavy elements is not fully reached even with 
half a million STAs. 

In the present work, we derive a formula for the spectrum of Configurationally-Resolved- 
Super-Ttransition-Arrays (CRSTA) which represents the exact sum of the spectra of all 
UTAs constituting the STA and sharing the same SCF solution. Out of the three assump- 
tions ((a)-(c) mentioned above), required by the traditional STA method, only the basic 
superconfiguration assumption (a) was used in the derivation of our new CRSTA method. 
As a consistency check we show that the radiation intensity, average energy and variance 
of the standard STAs are recovered from the CRSTA by an approximation based on a cu- 
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mulant expansion, truncated at the third term, i.e., the spectrum of a STA evaluated by 
the traditional method is the coarse-grained Gaussian approximation of the spectrum of the 
corresponding CRSTA. 

The plan of the manuscript is as follows: In section HIl some well known formulae required 
for the evaluation of the absorption coefficient in terms of the two-time dipole autocorre- 



lation functionj26|,[l2|],[27|,j28|] expanded in the eigenstates of the atomic Hamiltonian are 
summarized. The STA order of summation is briefly reviewed in section IIIII In section IIVI 
the formula for the CRSTA spectrum is derived. This formula is limited to the simple case 
where the Boltzmann factor is evaluated only with zero-order energies and the widths of 
the UTA are neglected (as in the original STA paper The derivation of the formula for 
the general case with the effect of electron-electron interaction in the Boltzmann factor and 
with the inclusion of the width of the UTA is derived in the Appendix. Section |V] contains 
a brief description of our new code, for the evaluation of the spectral absorption coefficient 
by both the standard STA and the new CRSTA methods, and a numerical example demon- 
strating the simplicity efficiency and accuracy of the new CRSTA method. A Summary and 
discussion are presented in section |Vll 

II. THE ABSORPTION COEFFICIENT 

The absorption coefficient, in hot dense plasmas in LTE, may be written in terms of the 



two-time autocorrelation function of the atomic many-electron dipole |12| . |26| . |27l| . |28j . |29| : 



/.at (E) = E'^^^n, Re {/,°° Ck (r) e^^^/^dr} . (1) 

/iat is defined as the fraction of the net absorbed radiation energy at energy E per unit 
radiation propagation length, uq is the atoms number density, Ck (t) is the two-time auto- 



correlation function of the atomic many-electron dipole 2^'■ 



Ck (r) = I3trace (^p^^d (r) d . (2) 

Eq. ([T]) is just one of many manifestations of the fluctuation-dissipation theorem connect- 
ing between the response of a given system to an external disturbance and the correlation 
of internal fluctuations of the system in the absence of the disturbance In Eq. ([2]), peq 
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is the equilibrium density matrix; d (r) is the the projection of the Heisenberg representa- 
tion of the atomic dipole operator, ^rj, on the polarization vector of the radiation field, 

i 

(averaging over all possible polarizations is implied), fi is the position operator of the i'th 
electron, d denotes the Kubo transform of d: 

d{(5) = ^ d{ihX)dX. (3) 

(For the relation between the correlation function Ck and more commonly used functions 
such as the symrnetrized correlation function C (r) = jgtrace {peq {d (r) d{0) + d (0) d{T))) , 



see references 



26| 



mm )• 

Expanding the dipole and the density operators in eigenfunctions of the atomic Hamil- 
tonian, H, with energies {En} using: 



^^e-/3(i3„-QM) ' (4) 



n 



where Q is the number of electrons, /i is the chemical potential, (3 = l/fc^T, T is the 
temperature and /c^ is the Boltzmann constant, the formula for /Xat becomes: 



/xat (E) = E^'finol'-^ Re d (r, /3) e-^^^l'W, (5) 

with 

^ (^, 13) = ye-P(E,-Q^^) 1^ |2 r[E,-E,y/k 

Y^e-f^i^^-Q^)^ ^ (6) 



and dif = {^\d\fj. Eq. (|5]) is equivalent to the Fermi golden rule in the form used in |l4||15l|. 
This can be easily seen by performing the r integral and using the relation 5 {Ef — Ei — E) = 
- lim^_>n Iin s -5 — „ . \. However, as will become clear from the results below, the 

TT 7 \ Ef—Ei—E—t-y j ' ' 

summation over states becomes an easier task when performed prior to the r integration. 



III. SUMMATION OVER CONFIGURATIONS 

Eq. (j6]) is a formula for ^ in terms of the exact many-electron energies and dipole matrix 
elements. As described in the introduction, in practice, these quantities are evaluated by a 
perturbation expansion, using the radial average potential approximation as the zero order. 
A configuration is a zero order many-electron state described by the occupation numbers 
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of the shells. Symbolically, a configuration is written as C = Y\ {^shisY' where a shell is 

s 

defined by the principal quantum number n^, the orbital angular momentum of the large 
component in the Dirac wave function, /g, and the total orbital+spin angular momentum 
is ; qs is the occupation number of the shell (configurations are degenerate states). Mutual 
electron-electron interaction lifts the degeneracy and splits the configurational energy. This 
effect is evaluated as a first order correction, i.e. as the sum of expectation values, in the zero 
order states of energies due to coulomb interaction between all pairs of electrons. In principle, 
to evaluate the first-order correction one should diagonalize the perturbation within each 
degenerate subspace. The methods discussed here involve only configurational averaged 
quantities, i.e. taking the trace, a process which does not require the diagonalization. 



Based on this picture, the STA method IJ] 15| splits the summation in Eq. ([6]) into 



a few stages. First, the spectrum is split into the different contributions of one-electron 
transitions, i.e.: 



a,b 

where i)"''' is the spectrum due to all possible transitions in which an electron transits from 
the shell nJaja to another shell nbhjb ■ The summation is over all configurations C in which 
the shell a has at least one electron i.e. the occupation number is > and the shell b 
has at least one hole, i.e. qt < 2jh + 1 . The summation over all accessible configurations is 
further partitioned by introducing an intermediate summation step over superconfigurations 



S(l|(l5|: 



(8) 



where N = ^^^^exp {Ec - Qfi)) = '^gcexp {-(3 (Ec - Qfx)). 

H cgh auc 
Each super configuration represents a particular distribution of the electrons between 

supershells ( a group of energetically adjacent atomic shells). 

In Eq.®, represents the contributions from the transition array C"*'', composed of 

all transitions from the shell nJaja in the configuration C to the shell rihlhjh, and gc is the 

zeroth order degeneracy. Following the STA method, we adopt the assumption of the UTA 

approach 1^ , that these transitions merge into an unresolved spectrum of a Gaussian shape 



(i.e a UTA). Thus, the three lowest energy moments of the UTA /g'' E^'' {AE^^) are used 
to construct the spectrum. In the context of the present work this means 

(r, P) = fc exp {-i {AE^')\' + tE^'r] (9) 



and: 



= — -J^ac^^p (Ec - Qfi)) rS (r, (3) 



J2fc9cexp{<!>-^{{qc},^,T',ZT)] 



(10) 



where 



^h' ({g^} , /3, r^ tr) ^ -(3 (Ec - Q/i) - | (AE^^) ' + ^E^V. (H) 

The moments of the UTA are represented in the following compact formulae [l^; The 
strength of a transition is: 

fc=<l^{9^-Q^)i{<^\\r\\b))\ (12) 

where (a ||r|| 6) is the reduced matrix element of the dipole. The configurational average 
of the energy is: 

Ec = + £1]^^ " ^'-'^ ^) ' (13) 

s s r 

with 

{s)^es + {s\-V{r)-f\s) (14) 

and 

(.,r) = FO(.,r)-i^^(l-5.Ao) ( ^ ^\ | GW(.,r), (15) 

\i/2 -1/2 y 

where F^^\G^^'^ are the Slater integrals corresponding to direct and exchange interaction, 

Js ^ Jr \ 

is the 3j Symbol, and Ig, k, l^. obey the triangle inequality. 

1/2 -1/2 j 

The center of gravity of the UTA is: 

E^' = D^' + 5^ (g. - 5sa) [of + ( - J^) ^ (a, h)) , (ig) 



where: 



= {b) - {a) , (17) 
Df = {{s,b) - {s,b)) (18) 



and 



^ia,b)^-j:g^g,{ ^ ^ ^ '\ I F(^) (a, 6) 



„¥i I 1 I V I -I / V I -i 

2 



even 



(19) 



The variance of the UTA is: 

(AEg^) ' = 5^ (g. - 5.a) i9s -qs- Ssb) {A'):' , (20) 

s 

where (A^)^^ is independentof the occupation numbers. 

These are all the building blocks necessary for the summation in Eq. (1T0|) . 

IV. CONFIGURATIONALLY-RESOLVED-SUPER-TRANSITION- ARRAYS 

For the simplicity of presentation we focus on the case in which the width of the UTA, 
(AE'g?)^, as well as the electron-electron interaction terms in the Boltzmann factor, i.e. in 
Ec (but not in E^^), are ignored. The treatment of the general case is deferred to the 
Appendix. Ignoring the width of the UTA, as well as the electron-electron interaction terms 
in Ec, Eq. (|TOl) is reduced to: 

^ab ^ 1 [S{Ec'Qt^) fab iE^b^/h 

^ E^ce-''(-C-«.)^^^ ' (21) 

Ces 

9s 



where the zeroth order degeneracy is: gc = Y[ \ \ ■ Fourier transform yields the 
STA spectrum. The standard STA method [ij] may be obtained by a cumulant expansion 



12| ) of Eq. (1211) and truncation at the third cumulant. This may be seen by writing the 



Taylor series for the factor e^^c''^/^''^ to obtain the expansion: 

oo 

^^ = mfiEh^'^/^r^'n\ (22) 



n=0 
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where: 



and fi = ^c,ce"^^-^^"^^)/c 



(23) 



Using the Gaussianity assumption, 



oo 
oo 



1 



/27r( Ae2 



exp 



yields the total radiation intensity, average energy and variance of the STA. Explicitly, 



(24) 



jab 



1 eab 



d 



di-ir/h) 
92 



In 



r=0 



i9g.''(r,/3) 



T = 



Mo 



T ill ?S5 



T=0 



/f2 

Moi — I 



Mr 

Mo' 



(25) 
(26) 
(27) 



_a(-iT/nY V -fS 

These results coincide with Eqs. (20), (21), (22) in reference [l^. In the practical application 
of the standard STA method, one is enforced to represent the spectrum by a large number 
of narrow Gaussian STAs, in order to minimize the error due to the truncation of the series. 
Typically, a STA run uses only a few self-consistent potentials (and sets of one-particle states 
and energies) for every degree of ionization, but a multitude of STAs. Thus, a multitude of 
STAs (and many more UTAs) share the same set of one-particle states and energies. 

Our new CRSTA method avoids the approximation of Gaussian STAs (i.e. the approx- 
imation in the truncation at the second cumulant). This is done by the application of the 
mathematical machinery of partition functions of independent particles directly to the con- 
tribution of transition arrays (actually their Fourier transform), to obtain the exact sum 
of all UTAs sharing the same one-particle states and energies. Explicitly, this is done by 
absorbing the time dependent exponent in the Boltzmann factor and using Eqs. (fT3l) . ( fT6|l 
to write the exponential factors in Eq. (J2T]) as: 



g^(,-KEc-Qf^)(,iE^''T/h 



n 



sec \ q. 



c 



) AJJ 



sec \ Q. 



c 



-gf{/3(£,-/.)+iDfr/n} 



(28) 
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where: 



(29) 



Next, we define the " superconfigurational degeneracy vector" with components (g)^ = Qs 
and the supertransitional "Complex Pseudo Partition Function" (CPPF): 



ceHsGC \ Qs 



and 



Using the well known combinatorial relations q 
















1 ^ H 



























^ 9 / 



(30) 



(31) 



and 



((a||.||b))^e-'«-"g'')-^'' 



9s — ^sa — ^sb 



sec 



q? - 



{Xf) 



((a||r||b))^e-«-"g'')-/'' 



(32) 



In Eq.f l32|) . the length of the vector 5a equals to the number of shells, where all components 
vanish except for the a component which has the value of 1. 

Note that the algebraic dependence of the formula for the partition function, Eq. ( !30l) . on 
the degeneracies, shell occupation numbers and X is the same as of the standard partition 
function. Therefore, it obey s the same recursion relations and is accessible to the efficient 
evaluation methods 

H,[l9l,Q,|2ll- Unlike references QQ 



19| 



20| 



where the 



partition function is used for the evaluation of the energy-moments of the spectrum, Eqs. 
( 132|) when Fourier Transformed with respect to time expresses the STA spectrum itself as a 
partition function. 
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Finally, the formula for the spectral absorption coefficient is obtained by using Eqs. ( 132|) 
and dSD: 

S,a,6 (33) 

V. NUMERICAL EXAMPLE 

We have written, from scratch, a new numerical code called CRSTA. This code can cal- 
elate the special abso.pt.o b. two optional methods, the standard STaH.Q .etKod 

and our new CRSTA method. In the STA method, single-particle energies are used to con- 
struct the real partition functions, which in turn are used together with the Slater integrals 
to construct the total intensity, average energy and variance of the STA. The spectrum 
of each STA is constructed as a Voigt function, which accounts for STA width as well as 
Doppler and electron impact effects. The absorption spectrum is obtained by the summa- 
tion of contributions from all STA's. In the CRSTA method, single-particle energies and 
Slater integrals are used to construct the CPPF at different times (Eq. ( 130|1 ). from which 
the spectral absorption coefficient is evaluated by Eq.f l33|) . 

The focus of the the ilustrating example was on the replacement of the traditional STA 
method of summation over configurations, which is based on the Gaussian approximation, 
by the CRSTA method in which the summation is exact (within the framework of the basic 
superconfiguration assumption (a)). Naturally, important technical and physical issues, 
which are common to the traditional and the new methods (such as stable and effecient 
evaluation of the partition function, and modeling of plasma effects ( 19|),(j30|)), where out 



of the focus of the present work. These issues where treated, in the present work, by the 
simplest possible methods. 

In both the new and traditinal methods, the relativistic single-particle radial eigenfunc- 
tions and eigenvalues are obtained from a self consistent solution (SCF) of the radial Dirac 
equation with the Hartree-Fock-Slater (HFS) potential in an ion-sphere. In our illustrating 
example, we use a SCF procedure, based on a very simple plasma model, that is briefly 
decribed as follow. In the flrst stage, an AA SCF calculation is done, providing the chem- 
ical potential and the average ionized electrons number. The chemical potential enforces 
neutrality within the AA's ion-sphere cell, where the free electrons density is approximated 
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by Thomas- Fermi model. In the second stage, a single SCF is calculated for each Q state 
(containing a single supershell and superconfiguration). For simplicity, we model each Q 
state as an ion-sphere that contains Q electrons and embedded in a uniform free electrons 
density, po- The former is calculated as the AA's ionized electrons number divided by the 
AA's cell volume. The neutrality of each Q state's ion-sphere is obtained by choosing its 
radius, rg, as: Q + 4/3po7rrQ = Z. 

The STA method requires two refinement loops. In the external refinement loop, the 
number of SCs sharing one-particle solutions with the same HFS potential is decreased. For 
example, in the AA limit, only one potential is used for all possible configurations sharing 
the same potential. In reference 1J| a different potential is used for each degree of ionization. 
Further refinement and convergence is obtained when few potentials are used for each degree 
of ionization. Yet, further spectral resolution is obtained by an additional internal loop, in 
which the number of Gaussian STAs sharing the same HFS solution is increased. In the 
CRSTA method, the internal loop of refinement is not necessary, since the exact analytical 
sum of contributions from all configurations sharing the same potential is represented by a 
single CPPF. On the other hand, the CRSTA method requires evaluation of CPPF on a time 
grid. We calulate the CPPFs using the fast recursion formulas of Bar-Shalom et al. [14] , by 
substituting the complex quantity Xf' {f3, r) rather than the standard Xs (= X^^ 0)). 
Of course, the complex factor of X"^ r) does not affect the numerical stabilty / instability 
of these formulas, therefore, one can use the same alternative stable methods 191]- 22], by 
substituting Xg X^^(/3,r), when numerical instabilities are expected (in our numerical 
example, we freeze several lowest energy shells to be fully occupied, and account small 
electrons number to be active in several supershells, therefore, insability is not expected). 
In our calculation, for each SCF potential, Vj, and allowed transition, aj — > bj, where 
ttj and bj are bound ionic shells we define a time grid. The resolution and size of the 
time grid are determined from the zeroth order transition energies , and an estimation of 
the expected width, F. Rapid oscillations in the integrand are elliminated by utilizing 
the fact that the Fourier transform (FT) in Eq. ( l33l) equals to the FT of the function 
Xf{T)U§^Q^^ (g-6a- Sb, (3, , shifted by the energy D^^ - Df. The time resolution. At, 
is estimated as 27r/At = n x — £:q,| x w, and the time interval as [0, 2m/r], where typically 
2<n<5, m>A and 0.1<w<0.3. In order to incorporate impact broadening effects, 
we also multiply by the factor exp(— r/2r) prior to the Fourier Transform. Typically, only 
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some 100-1000 CPPF calculations are required for each transition aj — )■ hj^m. order to obtain 
a UTA like resolution of the spectrum. 

Figure (1) shows the b-b absorption spectrum of the 8^3/2 — 4^5/2 transition in a Gold 
plasma at temperature of IkeV and density of 0.5(7r/cm^,(see similar example in p. 227 in 
Ref. ItI). Only zero order energies are used in the evaluation of the partition function and 

n r 

the widths of the UTAs are ignored [lj| (See however reference [15| and Appendix IVIip . 
Also, for the sake of simplicity we have used a single potential for each degree of ionization, 
having the ionic states Q = 11. .21 (shells below 3s shell remain fully occupied for these 
Q-states). The time grids contain 2000 grid points and resolution of O.lSa.u. The electron 
impact parameter, F ^ 0.04a.M., is calculated from the simple model of section 7.1.2 in Ref. 
Q]. For comparison, we have also evaluated the spectrum by the STA method with increasing 
number of supershells, and correspondingly increasing number of superconfigurations and 
Gaussian STAs. Figure (2) focuses on the energy range 3420-3460 eV. Clearly, a convergence 
is obtained only when a huge number of STAs (half a million) is used. Note that in the range 
3420-3430 eV even half a million of STAs do not fully reveal the spectral details. In this case, 
half a million partition functions and half a million recursion formulas were calculated for the 
energy and variance of each STA. In addition, half a million Voigt profiles were calculated 
for each STA. On the other hand, in the CRSTA calculation, only 22, 000 CPPF and 11 FTs 
were calculated. In order to check the analytical convergence of the STA spectrum to the 
CRSTA spectru, we have also calculated an extreme case of a UTA spectrum of a low ionic 
state, Q = 14, where only four electrons out of the 14 remain unfrozen and play an active 
combinatorial role. For this check, the STA spectrum is calculated as a sum of Lorentzian 
profiles, rather than Gaussians. The CRSTA and STA spectral profiles in this extreme case 
(not shown here) completely coincide. 
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Figure(l): b-b absorption spertrum of gold at T = IkeV and p = 0.5 g/cmP 
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Figure(2) Zoom of the energy range 3420-3460 eV in Figure(l) 



VI. SUMMARY AND DISCUSSION 

The main resuh of the present work is in Eqs. fl32|) and fl36l) . fl37|) for the evaluation of the 
spectrum of a STA in terms of a single CPPF. In order to analyze the difference between 
the new CRSTA method and the traditional STA method let us examine the derivation 
starting from Eq. (l^Tj) . This equation contains an instruction to sum over all the transitions 
constituting the STA, S, to get the Fourier transform of the STA spectrum , 'd^ij^P). 
The traditional STA formula for d'^ is obtained by Taylor expansion with respect to r, (Eq. 

15 



(122|) ). The Gaussian assumption (Eq. (12^ ). leads to the truncations of the expansion at the 
third term yielding formulae fl25l) . fl26|) and fl27j) for the radiation intensity, average energy 
and variance from which the STA spectrum is constructed as a Gaussian. The truncation 
of the Taylor expansion at the third term is justified in the short time limit. In the energy 
domain, it means a coarse graining of detailed structures finer than the variance of the 
STA. In contrast, in the CRSTA method, the summation to obtain d'^ (r, /3) is performed 
directly by the combinatorial steps fl28l) - fl32l) . which yield the exact result in terms of a single 
CPPF (Eqs. (!32|) and in the general case equation (l36l) . (l3711 ). Out of the three assumptions, 
required for the derivation of the traditional STA method [The basic superconfiguration 
assumption (a), the high-temperature approximation (b) and the unresolved supertransition 
array assumption (c)), only the first one is used for the derivation of the CRSTA method, 
i.e. this new method utilizes the simplicity of the analytical manipulations, enabled by the 
basic superconfiguration assumption (a), without suffering from the spectral coarse graining 
imposed by the unresolved supertransition array assumption (c). This is the reason why 
in the numerical example, presented in figures (l)-(2), a detailed spectrum, which required 
half a million narrow STAs for its resolution, is resolved by a few CRSTAs. 

The traditional STA concept was originally developed for plasmas at LTE. Later on, the 



idea was adapted also to the treatment of non-LTE plasma conditions (e.g. [32],|33|), and 
also to the treatment of the electronic degrees of freedom in the equation of state (e.g. 34|). 
The CRSTA method may be adapted also to these tasks. 

The possibility to extend the method to resolutions beyond the UTA should also be 
explored. Another direction which should be explored is the incorporation of the CRSTA 



method with screened hydrogenic model (SHM e.g.j35|) into a code for rapid (possibly in- 
line) evaluation of opacity. The CRSTA method removes the calculational bottleneck of 
summation over transitions, while the SHM removes the calculational bottleneck of the SCF 
process and the evaluation of Slater integrals. 
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VII. APPENDIX A: A FORMULA FOR THE SPECTRAL ABSORPTION COEF- 
FICIENT ACCOUNTING FOR UTA WIDTH AND FIRST ORDER CORRECTION 
IN THE BOLTZMANN FACTOR. 

The occupation numbers may be written as a vector of length of the number of shells, 

-^sheli: 

Using this definition and the explicit form of the configurational average energy, center of 
gravity and variance of a UTA (Eqs. f|T3|) . f|T6|) and f l20|) respectively), the scalar exponent 
^'c } 5 /^) ''"^i '^'^) (Eq- ( [TT]) ) may be written as a sum of three scalars; a scalar 0'^'' which 
is independent of g^, a scalar product between (f^ and a vector of coefficients T"^ which is 
independent of g^, and a quadratic form in the vector (f^ with a real symmetric iVsheii * ^sheii 
matrix of coefficients, Q which is independent of g^. 

The quadratic form prevents a direct application of the combinatorial manipulations of 

^ y 

Eqs ( 128|) - (!32|) . To cure this problem we diagonalize the matrix of coefficients, Q and 



apply the Hubbard- Stratonovich transformation 



24| (Eq. flM|) below) which eliminates the 



nonlinear dependence on occupation numbers. Explicitly, this is done as follows: 
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First, the quadratic form is evaluated in a rotated system; 



which is chosen so that the matrix Q is diagonal. Explicitly, the elements of the k 



eigencector of Q obey: 



i.e. Q is diagonalized by a matrix with the elements R. 



rk ; 



A I — Rjr^rsRsk Rjr^kRrk 
jk g^j. ^ 

= ^kYl Rjr-Rrk = ^k^jk, 



and Q is obtained from A by the inverse of this symmetric transformation: 



Using Eq. ( 1121) for the dipole matrix element together with the integral identity: 

exp (p?)' Af ) = ^ exp (^-1 |c,| + ixp^' ^2Xf \cr\^ dx 
the sum in Eq. (|TOl) is written as: 



ab 



with 



= E {fc9c exp |-0 - E gf T,} } A^^' 

= n {a/I^ X^L^Xp \Cr\xl + iXrV^\/2\f \Cr\^ dx^ 

I-oo - Too (--2 E \Cr\ X? + ^ E ?f E XrR:t./2Xf^\) n y^dXr 

\ r s r J T 



V t s s 

The matrices ^ and P and the vector ,^in Eq.f l35p are defined by 

fr^) =Rf,./2XfWl 

\ / sr 



(34) 



(35) 
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and: 



Next, define: 



St 



i = T -x. 



-1 



rt 



T^' {i) = exp (-i E E PstU^ {^) ' 



1 1 V 27r • 

r 

With this definition, the summation over configuration takes the form: 



X 



= E { exp I - - E T, 1 1 

n - n exp (-1 E E + ^ E e.) I C^^) " 

= exp {—6} {{a I |r| I 6))^ 



(36) 



see \ q. 



where the expression gc = Y[ I ^'^ I ^"^^ configurational degeneracy was used. 



Now define 

Xf {13, t) = exp {-Tf + z^f } 

and the complex pseudo partition function: 



fj-ab 



En, 
E n 



(37) 



9^ \ f^ab^'^"^ 



C{Sum(gC')=Q}seC ^ Qs 

and apply the same combinatorial steps as in Eq. (132|) . The result is a the generalization of 
Eq. m- 
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/^at (E) = f0no (1 - e-f^^) i{a\\r\\b)fg^g, 

X Re dr {exp { -6^^ (r) } rf^f {^f (f) f/i^g-i (^^ " I " 5., r, tj (t) } } • 

(38) 

Eq.( l38|) is a formula for the spectral absorption coefficient which accounts for the width of 
the UTA, as well as the first order correction in the Boltzmann factor without the restriction 
of the high temperature approximation. 
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